{ "cells": [ { "cell_type": "markdown", "metadata": { "id": "N2IUDwhPvDD2" }, "source": [ "# 1. Prepare input data of Kano River\n", "*Written by Men Vuthy, 2022*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Import modules**" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "executionInfo": { "elapsed": 3402, "status": "ok", "timestamp": 1647337073216, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "GUtVEQQsYOqG" }, "outputs": [], "source": [ "import os\n", "import numpy as np\n", "import pandas as pd\n", "\n", "import matplotlib.pyplot as plt\n", "\n", "import rasterio\n", "import geopandas as gpd" ] }, { "cell_type": "markdown", "metadata": { "id": "4kzODWJeteuT" }, "source": [ "**Label**" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "executionInfo": { "elapsed": 442, "status": "ok", "timestamp": 1647337078596, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "0eRWOvq2YlB1" }, "outputs": [], "source": [ "kano_class_img = rasterio.open('data/kano_river/out_img/class/kano_class.tiff')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "executionInfo": { "elapsed": 762, "status": "ok", "timestamp": 1647337079355, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "8WuGITVVfVpR" }, "outputs": [], "source": [ "kano_class = kano_class_img.read(1).reshape(-1)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "executionInfo": { "elapsed": 3, "status": "ok", "timestamp": 1647337079356, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "34uDVmPonYvn" }, "outputs": [], "source": [ "kano_label = pd.DataFrame({'label':kano_class})" ] }, { "cell_type": "markdown", "metadata": { "id": "efyxwD0c8I8z" }, "source": [ "**Dataframe of Season 1**" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "executionInfo": { "elapsed": 10, "status": "ok", "timestamp": 1647337079750, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "dxzPDKy_fvWT" }, "outputs": [], "source": [ "season1_rgbn = rasterio.open('data/kano_river/out_img/rgbn/kano_20180429_rgbn.tiff')\n", "season1_ndvi = rasterio.open('data/kano_river/out_img/ndvi/kano_20180429_ndvi.tiff')\n", "season1_ndwi = rasterio.open('data/kano_river/out_img/ndwi/kano_20180429_ndwi.tiff')\n", "season1_bsi = rasterio.open('data/kano_river/out_img/bsi/kano_20180429_bsi.tiff')" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "executionInfo": { "elapsed": 2865, "status": "ok", "timestamp": 1647337082606, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "A9aRPiqjrBjC" }, "outputs": [], "source": [ "blue_1 = season1_rgbn.read(1).reshape(-1)\n", "green_1 = season1_rgbn.read(2).reshape(-1)\n", "red_1 = season1_rgbn.read(3).reshape(-1)\n", "nir_1 = season1_rgbn.read(4).reshape(-1)\n", "\n", "ndvi_1 = season1_ndvi.read(1).reshape(-1)\n", "ndwi_1 = season1_ndwi.read(1).reshape(-1)\n", "bsi_1 = season1_bsi.read(1).reshape(-1)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "executionInfo": { "elapsed": 696, "status": "ok", "timestamp": 1647337084966, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "cKMuZhn8rJgR" }, "outputs": [], "source": [ "Season_1 = pd.DataFrame({'B1':blue_1, 'G1':green_1, 'R1':red_1, 'NIR1':nir_1, 'NDVI1':ndvi_1, 'NDWI1':ndwi_1, 'BSI1':bsi_1})" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 423 }, "executionInfo": { "elapsed": 5, "status": "ok", "timestamp": 1647337084966, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "y8dt2IOysWrh", "outputId": "15306d04-392f-4580-e0a9-738703332b9c" }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
B1G1R1NIR1NDVI1NDWI1BSI1
00000NaNNaNNaN
10000NaNNaNNaN
20000NaNNaNNaN
30000NaNNaNNaN
40000NaNNaNNaN
........................
236637960000NaNNaNNaN
236637970000NaNNaNNaN
236637980000NaNNaNNaN
236637990000NaNNaNNaN
236638000000NaNNaNNaN
\n", "

23663801 rows × 7 columns

\n", "
" ], "text/plain": [ " B1 G1 R1 NIR1 NDVI1 NDWI1 BSI1\n", "0 0 0 0 0 NaN NaN NaN\n", "1 0 0 0 0 NaN NaN NaN\n", "2 0 0 0 0 NaN NaN NaN\n", "3 0 0 0 0 NaN NaN NaN\n", "4 0 0 0 0 NaN NaN NaN\n", "... .. .. .. ... ... ... ...\n", "23663796 0 0 0 0 NaN NaN NaN\n", "23663797 0 0 0 0 NaN NaN NaN\n", "23663798 0 0 0 0 NaN NaN NaN\n", "23663799 0 0 0 0 NaN NaN NaN\n", "23663800 0 0 0 0 NaN NaN NaN\n", "\n", "[23663801 rows x 7 columns]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Season_1" ] }, { "cell_type": "markdown", "metadata": { "id": "5tIbL4BkpJa2" }, "source": [ "**Dataframe of Season 2**" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "executionInfo": { "elapsed": 4, "status": "ok", "timestamp": 1647337084966, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "yYKDT-cLpJJr" }, "outputs": [], "source": [ "season2_rgbn = rasterio.open('data/kano_river/out_img/rgbn/kano_20180715_rgbn.tiff')\n", "season2_ndvi = rasterio.open('data/kano_river/out_img/ndvi/kano_20180715_ndvi.tiff')\n", "season2_ndwi = rasterio.open('data/kano_river/out_img/ndwi/kano_20180715_ndwi.tiff')\n", "season2_bsi = rasterio.open('data/kano_river/out_img/bsi/kano_20180715_bsi.tiff')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "executionInfo": { "elapsed": 2410, "status": "ok", "timestamp": 1647337089748, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "lLq3mThwpZnI" }, "outputs": [], "source": [ "blue_2 = season2_rgbn.read(1).reshape(-1)\n", "green_2 = season2_rgbn.read(2).reshape(-1)\n", "red_2 = season2_rgbn.read(3).reshape(-1)\n", "nir_2 = season2_rgbn.read(4).reshape(-1)\n", "\n", "ndvi_2 = season2_ndvi.read(1).reshape(-1)\n", "ndwi_2 = season2_ndwi.read(1).reshape(-1)\n", "bsi_2 = season2_bsi.read(1).reshape(-1)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "executionInfo": { "elapsed": 420, "status": "ok", "timestamp": 1647337090159, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "RAqkedgSpkXn" }, "outputs": [], "source": [ "Season_2 = pd.DataFrame({'B2':blue_2, 'G2':green_2, 'R2':red_2, 'NIR2':nir_2, 'NDVI2':ndvi_2, 'NDWI2':ndwi_2, 'BSI2':bsi_2})" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 423 }, "executionInfo": { "elapsed": 9, "status": "ok", "timestamp": 1647337090160, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "j2Zp2QI2pru4", "outputId": "1860ce19-dee3-4c61-9666-542ca06f43ae" }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
B2G2R2NIR2NDVI2NDWI2BSI2
00000NaNNaNNaN
10000NaNNaNNaN
20000NaNNaNNaN
30000NaNNaNNaN
40000NaNNaNNaN
........................
236637960000NaNNaNNaN
236637970000NaNNaNNaN
236637980000NaNNaNNaN
236637990000NaNNaNNaN
236638000000NaNNaNNaN
\n", "

23663801 rows × 7 columns

\n", "
" ], "text/plain": [ " B2 G2 R2 NIR2 NDVI2 NDWI2 BSI2\n", "0 0 0 0 0 NaN NaN NaN\n", "1 0 0 0 0 NaN NaN NaN\n", "2 0 0 0 0 NaN NaN NaN\n", "3 0 0 0 0 NaN NaN NaN\n", "4 0 0 0 0 NaN NaN NaN\n", "... .. .. .. ... ... ... ...\n", "23663796 0 0 0 0 NaN NaN NaN\n", "23663797 0 0 0 0 NaN NaN NaN\n", "23663798 0 0 0 0 NaN NaN NaN\n", "23663799 0 0 0 0 NaN NaN NaN\n", "23663800 0 0 0 0 NaN NaN NaN\n", "\n", "[23663801 rows x 7 columns]" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Season_2" ] }, { "cell_type": "markdown", "metadata": { "id": "k_HvDVLoqCGI" }, "source": [ "**Dataframe of Season 3**" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "executionInfo": { "elapsed": 463, "status": "ok", "timestamp": 1647337092434, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "_x0caSaXptHI" }, "outputs": [], "source": [ "season3_rgbn = rasterio.open('data/kano_river/out_img/rgbn/kano_20181029_rgbn.tiff')\n", "season3_ndvi = rasterio.open('data/kano_river/out_img/ndvi/kano_20181029_ndvi.tiff')\n", "season3_ndwi = rasterio.open('data/kano_river/out_img/ndwi/kano_20181029_ndwi.tiff')\n", "season3_bsi = rasterio.open('data/kano_river/out_img/bsi/kano_20181029_bsi.tiff')" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "executionInfo": { "elapsed": 4095, "status": "ok", "timestamp": 1647337097114, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "7rdZVc-Mqige" }, "outputs": [], "source": [ "blue_3 = season3_rgbn.read(1).reshape(-1)\n", "green_3 = season3_rgbn.read(2).reshape(-1)\n", "red_3 = season3_rgbn.read(3).reshape(-1)\n", "nir_3 = season3_rgbn.read(4).reshape(-1)\n", "\n", "ndvi_3 = season3_ndvi.read(1).reshape(-1)\n", "ndwi_3 = season3_ndwi.read(1).reshape(-1)\n", "bsi_3 = season3_bsi.read(1).reshape(-1)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "executionInfo": { "elapsed": 576, "status": "ok", "timestamp": 1647337097686, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "cUdE7Gq3qpVX" }, "outputs": [], "source": [ "Season_3 = pd.DataFrame({'B3':blue_3, 'G3':green_3, 'R3':red_3, 'NIR3':nir_3, 'NDVI3':ndvi_3, 'NDWI3':ndwi_3, 'BSI3':bsi_3})" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 423 }, "executionInfo": { "elapsed": 5, "status": "ok", "timestamp": 1647337097687, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "ETLVsonaq0qo", "outputId": "06530c17-56f1-4892-c83e-c9f67d153c37" }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
B3G3R3NIR3NDVI3NDWI3BSI3
00000NaNNaNNaN
10000NaNNaNNaN
20000NaNNaNNaN
30000NaNNaNNaN
40000NaNNaNNaN
........................
236637960000NaNNaNNaN
236637970000NaNNaNNaN
236637980000NaNNaNNaN
236637990000NaNNaNNaN
236638000000NaNNaNNaN
\n", "

23663801 rows × 7 columns

\n", "
" ], "text/plain": [ " B3 G3 R3 NIR3 NDVI3 NDWI3 BSI3\n", "0 0 0 0 0 NaN NaN NaN\n", "1 0 0 0 0 NaN NaN NaN\n", "2 0 0 0 0 NaN NaN NaN\n", "3 0 0 0 0 NaN NaN NaN\n", "4 0 0 0 0 NaN NaN NaN\n", "... .. .. .. ... ... ... ...\n", "23663796 0 0 0 0 NaN NaN NaN\n", "23663797 0 0 0 0 NaN NaN NaN\n", "23663798 0 0 0 0 NaN NaN NaN\n", "23663799 0 0 0 0 NaN NaN NaN\n", "23663800 0 0 0 0 NaN NaN NaN\n", "\n", "[23663801 rows x 7 columns]" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Season_3" ] }, { "cell_type": "markdown", "metadata": { "id": "S2-5njKqq-Ou" }, "source": [ "**Dataframe of Season 4**" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "executionInfo": { "elapsed": 4, "status": "ok", "timestamp": 1647337097687, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "RiTAGhkjq3T4" }, "outputs": [], "source": [ "season4_rgbn = rasterio.open('data/kano_river/out_img/rgbn/kano_20190129_rgbn.tiff')\n", "season4_ndvi = rasterio.open('data/kano_river/out_img/ndvi/kano_20190129_ndvi.tiff')\n", "season4_ndwi = rasterio.open('data/kano_river/out_img/ndwi/kano_20190129_ndwi.tiff')\n", "season4_bsi = rasterio.open('data/kano_river/out_img/bsi/kano_20190129_bsi.tiff')" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "executionInfo": { "elapsed": 2876, "status": "ok", "timestamp": 1647337100559, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "7Tb5TVotrQc4" }, "outputs": [], "source": [ "blue_4 = season4_rgbn.read(1).reshape(-1)\n", "green_4 = season4_rgbn.read(2).reshape(-1)\n", "red_4 = season4_rgbn.read(3).reshape(-1)\n", "nir_4 = season4_rgbn.read(4).reshape(-1)\n", "\n", "ndvi_4 = season4_ndvi.read(1).reshape(-1)\n", "ndwi_4 = season4_ndwi.read(1).reshape(-1)\n", "bsi_4 = season4_bsi.read(1).reshape(-1)" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "executionInfo": { "elapsed": 8, "status": "ok", "timestamp": 1647337100560, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "Ly1uhXo4rXX3" }, "outputs": [], "source": [ "Season_4 = pd.DataFrame({'B4':blue_4, 'G4':green_4, 'R4':red_4, 'NIR4':nir_4, 'NDVI4':ndvi_4, 'NDWI4':ndwi_4, 'BSI4':bsi_4})" ] }, { "cell_type": "code", "execution_count": 20, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 423 }, "executionInfo": { "elapsed": 6, "status": "ok", "timestamp": 1647337100561, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "p_ivo4ZordUW", "outputId": "fcae1564-59e2-4ade-8e3f-04f551f688b2" }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
B4G4R4NIR4NDVI4NDWI4BSI4
00000NaNNaNNaN
10000NaNNaNNaN
20000NaNNaNNaN
30000NaNNaNNaN
40000NaNNaNNaN
........................
236637960000NaNNaNNaN
236637970000NaNNaNNaN
236637980000NaNNaNNaN
236637990000NaNNaNNaN
236638000000NaNNaNNaN
\n", "

23663801 rows × 7 columns

\n", "
" ], "text/plain": [ " B4 G4 R4 NIR4 NDVI4 NDWI4 BSI4\n", "0 0 0 0 0 NaN NaN NaN\n", "1 0 0 0 0 NaN NaN NaN\n", "2 0 0 0 0 NaN NaN NaN\n", "3 0 0 0 0 NaN NaN NaN\n", "4 0 0 0 0 NaN NaN NaN\n", "... .. .. .. ... ... ... ...\n", "23663796 0 0 0 0 NaN NaN NaN\n", "23663797 0 0 0 0 NaN NaN NaN\n", "23663798 0 0 0 0 NaN NaN NaN\n", "23663799 0 0 0 0 NaN NaN NaN\n", "23663800 0 0 0 0 NaN NaN NaN\n", "\n", "[23663801 rows x 7 columns]" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Season_4" ] }, { "cell_type": "markdown", "metadata": { "id": "LhqoBmmMvczS" }, "source": [ "**Features and Label**" ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "executionInfo": { "elapsed": 1824, "status": "ok", "timestamp": 1647337105978, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "-MMyjJRkrfjX" }, "outputs": [], "source": [ "# concat dataframe of all seasons and label into one dataframe\n", "Kano_dataframe = pd.concat([Season_1, Season_2, Season_3, Season_4, kano_label], axis=1)" ] }, { "cell_type": "code", "execution_count": 22, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 423 }, "executionInfo": { "elapsed": 5012, "status": "ok", "timestamp": 1647334593379, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "mgH7H4PluPw3", "outputId": "4915d489-1699-4e28-bf91-0541fde875c8" }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
B1G1R1NIR1NDVI1NDWI1BSI1B2G2R2...NDWI3BSI3B4G4R4NIR4NDVI4NDWI4BSI4label
00000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
10000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
20000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
30000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
40000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
..................................................................
236637960000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
236637970000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
236637980000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
236637990000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
236638000000NaNNaNNaN000...NaNNaN0000NaNNaNNaN0
\n", "

23663801 rows × 29 columns

\n", "
" ], "text/plain": [ " B1 G1 R1 NIR1 NDVI1 NDWI1 BSI1 B2 G2 R2 ... NDWI3 BSI3 \\\n", "0 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "1 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "2 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "3 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "4 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "... .. .. .. ... ... ... ... .. .. .. ... ... ... \n", "23663796 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "23663797 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "23663798 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "23663799 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "23663800 0 0 0 0 NaN NaN NaN 0 0 0 ... NaN NaN \n", "\n", " B4 G4 R4 NIR4 NDVI4 NDWI4 BSI4 label \n", "0 0 0 0 0 NaN NaN NaN 0 \n", "1 0 0 0 0 NaN NaN NaN 0 \n", "2 0 0 0 0 NaN NaN NaN 0 \n", "3 0 0 0 0 NaN NaN NaN 0 \n", "4 0 0 0 0 NaN NaN NaN 0 \n", "... .. .. .. ... ... ... ... ... \n", "23663796 0 0 0 0 NaN NaN NaN 0 \n", "23663797 0 0 0 0 NaN NaN NaN 0 \n", "23663798 0 0 0 0 NaN NaN NaN 0 \n", "23663799 0 0 0 0 NaN NaN NaN 0 \n", "23663800 0 0 0 0 NaN NaN NaN 0 \n", "\n", "[23663801 rows x 29 columns]" ] }, "execution_count": 22, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Kano_dataframe" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "executionInfo": { "elapsed": 965, "status": "ok", "timestamp": 1647337335648, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "z9jaUKUfv10G" }, "outputs": [], "source": [ "# Remove all rows where NDVI1 is NaN\n", "Kano_classified = Kano_dataframe.dropna(subset=['NDVI1'])" ] }, { "cell_type": "code", "execution_count": 24, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 423 }, "executionInfo": { "elapsed": 399, "status": "ok", "timestamp": 1647337338259, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "zF0sqf5XzBRE", "outputId": "7a4cc005-48a4-4c4c-dd89-e7ed1cefef8b" }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
B1G1R1NIR1NDVI1NDWI1BSI1B2G2R2...NDWI3BSI3B4G4R4NIR4NDVI4NDWI4BSI4label
58463689099299519650.242947-0.2141660.303773860912819...-0.3280720.32548265574576811480.207700-0.1722380.3074626
58463790099099020080.255461-0.2254270.306076882936865...-0.3041470.33212965272273510720.195912-0.1541650.3109036
58463890197295818620.235430-0.1981410.307286877915854...-0.2935470.3294326367097339830.155239-0.1205630.3127326
58463980188986315760.205786-0.1606420.297341846857799...-0.2902260.3193626196837058740.116680-0.0808910.3147626
58464084289389613010.093848-0.0638360.315164821811749...-0.2346300.3131645726426527390.072273-0.0281230.3072356
..................................................................
2200043079494799018670.221129-0.2119060.298590835877864...-0.2293480.3268204105245559020.247357-0.2253990.2877256
220004319261057111019110.177388-0.1701950.309133104710821071...-0.2004320.3220064475395929210.226740-0.2218940.3089926
220004329261057111019110.177388-0.1701950.309133104710821071...-0.2004320.3220064475395929210.226740-0.2218940.3089926
220004339781124121219190.136559-0.1422560.313926123512691296...-0.1661990.3182654905936508350.134188-0.1281850.3077326
220047518921018107018190.171261-0.1644870.309384889937920...-0.1957600.3139583835165457930.194771-0.1709390.2757296
\n", "

692246 rows × 29 columns

\n", "
" ], "text/plain": [ " B1 G1 R1 NIR1 NDVI1 NDWI1 BSI1 B2 G2 \\\n", "584636 890 992 995 1965 0.242947 -0.214166 0.303773 860 912 \n", "584637 900 990 990 2008 0.255461 -0.225427 0.306076 882 936 \n", "584638 901 972 958 1862 0.235430 -0.198141 0.307286 877 915 \n", "584639 801 889 863 1576 0.205786 -0.160642 0.297341 846 857 \n", "584640 842 893 896 1301 0.093848 -0.063836 0.315164 821 811 \n", "... ... ... ... ... ... ... ... ... ... \n", "22000430 794 947 990 1867 0.221129 -0.211906 0.298590 835 877 \n", "22000431 926 1057 1110 1911 0.177388 -0.170195 0.309133 1047 1082 \n", "22000432 926 1057 1110 1911 0.177388 -0.170195 0.309133 1047 1082 \n", "22000433 978 1124 1212 1919 0.136559 -0.142256 0.313926 1235 1269 \n", "22004751 892 1018 1070 1819 0.171261 -0.164487 0.309384 889 937 \n", "\n", " R2 ... NDWI3 BSI3 B4 G4 R4 NIR4 NDVI4 \\\n", "584636 819 ... -0.328072 0.325482 655 745 768 1148 0.207700 \n", "584637 865 ... -0.304147 0.332129 652 722 735 1072 0.195912 \n", "584638 854 ... -0.293547 0.329432 636 709 733 983 0.155239 \n", "584639 799 ... -0.290226 0.319362 619 683 705 874 0.116680 \n", "584640 749 ... -0.234630 0.313164 572 642 652 739 0.072273 \n", "... ... ... ... ... ... ... ... ... ... \n", "22000430 864 ... -0.229348 0.326820 410 524 555 902 0.247357 \n", "22000431 1071 ... -0.200432 0.322006 447 539 592 921 0.226740 \n", "22000432 1071 ... -0.200432 0.322006 447 539 592 921 0.226740 \n", "22000433 1296 ... -0.166199 0.318265 490 593 650 835 0.134188 \n", "22004751 920 ... -0.195760 0.313958 383 516 545 793 0.194771 \n", "\n", " NDWI4 BSI4 label \n", "584636 -0.172238 0.307462 6 \n", "584637 -0.154165 0.310903 6 \n", "584638 -0.120563 0.312732 6 \n", "584639 -0.080891 0.314762 6 \n", "584640 -0.028123 0.307235 6 \n", "... ... ... ... \n", "22000430 -0.225399 0.287725 6 \n", "22000431 -0.221894 0.308992 6 \n", "22000432 -0.221894 0.308992 6 \n", "22000433 -0.128185 0.307732 6 \n", "22004751 -0.170939 0.275729 6 \n", "\n", "[692246 rows x 29 columns]" ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Kano_classified" ] }, { "cell_type": "code", "execution_count": 25, "metadata": { "executionInfo": { "elapsed": 13034, "status": "ok", "timestamp": 1647337354951, "user": { "displayName": "MEN Vuthy", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Ggz_fgR6NAqa_UnWAIIDjX9pd1OkdeBn_u164cC=s64", "userId": "02183312642253897909" }, "user_tz": -540 }, "id": "kfDWqViT1pRI" }, "outputs": [], "source": [ "# Data dir\n", "data_dir = \"data/kano_river/out_img/classified\"\n", "\n", "# Output dataframe\n", "out_df = os.path.join(data_dir, 'kano_classified.csv')\n", "\n", "# export dataframe to csv\n", "Kano_classified.to_csv(out_df, index=False)" ] }, { "cell_type": "code", "execution_count": 26, "metadata": { "id": "ut1eM8opzDWU" }, "outputs": [], "source": [ "# Data dir\n", "data_dir = \"data/kano_river/out_img/classified\"\n", "\n", "# Output dataframe\n", "out_df = os.path.join(data_dir, 'kano_classified_index.csv')\n", "\n", "# export dataframe to csv\n", "Kano_classified.to_csv(out_df)" ] } ], "metadata": { "colab": { "authorship_tag": "ABX9TyPYU4+vk0i3H7YJwWPgbYAl", "name": "5. Input Data for Supervised-ML- Kano River.ipynb", "provenance": [] }, "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.10" } }, "nbformat": 4, "nbformat_minor": 4 }